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Abstract 



^ . This is a brief description of how to derive the local "atomic" potentials from the Self- 

Consistent Atomic Deformation (SCAD) model density function. Particular attention is paid 
. to the spherically averaged case. 



S : 1 The SCAD Model 

■ The Self-Consistent Atomic Deformation (SCAD) model has been described as an extension of 

the work of Gordon and Kim [^. In can also be viewed as an approximation to density functional 
theory, which is the point of view taken in this paper. 

The SCAD model approximates the exact electronic density functional 0] by the form: 



E[n{r)] =J2nhir)] + ^TtfHv)] -Y.TTF[ndr)]j + F[n{r)] . (1) 

Here: 

1. To[nj(r)] is the kinetic energy of non-interacting electrons centered about the site Rj, 

roK(r)] = - E / r^a{r)y'i^^air) d\ , (2) 

where the ipia are an orthonormal set of wave functions such that 

n.(r) = 5:i^,,(r)r (3) 

a 

(We use atomic units, with h = 2m = e^/2 = 1, where m is the mass of the electron and — e 
is its charge.) 
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2. We assume that the charge associated with each site, 

Ni= [ n,{r)Sr , (4) 



is fixed. That is, there is no charge transfer between sites during the iteration of the problem 
to self-consistency. 

3. The total density of the system is given by 

n(r) = ^ni(r + R,) . (5) 

i 

4. The Thomas-Fermi kinetic energy is given (in Rydberg units) by 

TtfHt)] = ^(37r2)l J nt(r)dV . (6) 

The sum over atomic Thomas-Fermi kinetic energies in (||) is to cancel this contribution to 
E[n] from the individual atoms, where the kinetic energy is given by (Q). The Thomas- Fermi 
terms in parenthesis in (|l|) thus represent the overlap kinetic energy, and vanish when the 
atoms are separated by large distances. 

5. The functional -F[n(r)] contains the remaining density functional terms ^, including the 
Coulomb interaction terms (electron-electron, electron- nucleus, and nucleus-nucleus), which 
is a simple functional of the total density (^), and the exchange-correlation term, which is 
local functional of the density within the Local Density Approximation (LDA). The functional 
derivative of F with respect to the total density is just the Kohn-Sham potential Q, 

, , , , 6F[n{r)] 
..s[n(r);r] = ^iA^. (7) 

This potential includes both the Coulomb and the exchange-correlation potentials. 

2 Finding the potential 

In this section we wish to find a way to determine the individual charge densities nj(r). We begin 
by noting that if rij is, as we assume, "v-representable" , then there is a one-to-one correspondence 
between the density rij and a potential Vi, such that the wave functions described in (|2|) and (|3|) 
satisfy the Schrodinger equation 

- V^ViQ:(r) + Vi{r)ipia{r) = eiatpia{r) , (8) 

This is equivalent to rewriting the SCAD density Functional (ffl) in the form 



E[Mr)}] =J2To[v,{r)] + i^TrpHr)] - ^ rTi.[n,(r)] | + F[n(r)] . 



(9) 



2 



Now consider changing the potential in (^) 

Viir) ^ v,ir) + 6viir) . (10) 
This change in the potential at the site i changes the associated wave functions 

-fpiair) ipiair) + 6lpia{r) , (11) 

where, to first order in 6vi, 

ha + - Vi{r)]Stpia{r) = {Svi{r) - 5eia)V'ia(r) . (12) 
Without loss of generality we may take each Sip to be orthogonal to the corresponding ip: 

5V*„(r) Via(r) = 0. (13) 



Multiplying both sides of (|T^ by "01*0: (r) and integrating, we obtain the first order change in the 
eigenvalue: 

6eia= f \i'ia{r)\^6vi{r)cfr . (14) 



The electron density on site i then changes by an amount 

5n,(r) = 25R5]#*,(r) V'ia(r) , (15) 

a 

which is a change of 0[5t;j]. (JR.z is the real part of z.) Note that the densities rij on the other sites 
will not change, since by supposition we are only changing the spherical potential on site i, and 
there is no charge transfer. Thus 

5n{r) = 6ni{r + -Ri) . (16) 
The kinetic energy Tq (^) changes by an amount 

6To = -2^Y. [ #L(r)V2v^i,(r) . (17) 

By this becomes 

(r)]^i„(r)(iV . (18) 

Applying ( |l^ ) and (15), we find 

6To = - [vi{r)6ni{r) . (19) 



The other parts of the density functional change in a straight-forward way: 

6F[n{r)] = [ VKs[n{r);r]6ni{r) d^r , (20) 
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and 



where 



STTpinir)] = J VTF[n{r);r] 6ni{r)(fr , 

6TTF[ni{r)] = J VTF[n{r); r] 5ni{r) d^r , 
STtfHt)] 



is the Thomas-Fermi potential. 



6n{r) 



(21) 
(22) 
(23) 



Substituting equations (19-23) into the energy formula (|^, we find that it changes by an amount 
5E[{vi{r)}] = / {vKs[n{r); r] + VTF[n{r); r] - VTF[ni{r); r] - w;(r)} 6ni{r) d^r . (24) 



Equation ( p^ must hold for an arbitrary infinitesimal change in the potential 5vi. By (15) and 
([l^, this produces a corresponding change in 5ni. Hence (24) must hold for arbitrary infinitesimal 
changes in drij, which can only occur if the term in curly brackets vanishes, i.e., if 



Vi{v) = VKs[n{r); r] + VTF[n{r); r] - VTF[ni{v)] r] 



(25) 



3 Spherical Averages 

When using the SCAD model for crystals with atoms at high symmetry sites, it is often useful to 
assume that the densities rii are spherically symmetric. This can be achieved by assuming that the 
Hi represent a closed shell system, and that each potential Vi is spherical. Under these assumptions, 
the formalism developed in Section |2| remains the same until equation (|2j) , except that we replace 
the arbitrary potential Vi{r) by the spherical potential Vi{r), and restrict the potential changes to 
the spherical 5vi{r). At that point we note that because of the assumption that rii represents a 
closed shell, (5nj(r) must also be spherical. Thus (24) can be written 

5E[{vi{r)}\ = j {vKs[n{Y); r] + VTF[n{v); r] - i-tf ["-i (r) ; r] - Vi{r)} 6ni{r) d^r . (26) 

This holds for arbitrary spherical potential changes 6vi{r), whence Vi{r) must obey 

Vi{r) = j {'"Ks[n{r)] r] + VTF[n{r); r] - fTFK(r); r]} df] , (27) 

where the integral is over the surface of a sphere of radius r. Note that this is just the spherical 
average of the full potential (p^): 

v^{r) = ^ j Vi{v)dn . (28) 
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4 Discussion 

The essential difference between the SCAD model and the more familiar Kohn-Sham LDA formu- 
lation is that the SCAD model approximates kinetic energy in low-density regions (outside atomic 
cores) with the aid of a local functional. While the total density in the SCAD model is expressed as 
a sum of localized densities, this should not be a major source of error, provided that the localized 
densities are given adequate flexibility. Solving the atomic calculations self consistently, with the 
potential for each site determined from the densities and potentials of the rest of the system, puts 
the atomic energy levels on a common scale. Just as in the Kohn-Sham treatment, the lowest 
energy levels are occupied to achieve the minimum total energy. 

The SCAD model offers some advantages over the Kohn-Sham LDA approach. It is generally 
more efficient, particularly for large systems, because the computational labor increases linearly 
with the number of atoms in the system. Moreover, it is ideally suited for application on a massively 
parallel computer. Finally, we note that the SCAD model has the correct behavior in the limit 
of infinitely separated atoms, whereas, each electron in the Kohn-Sham model is delocalized and 
remains so in this limit. 
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